Quantum fields near phantom-energy 'sudden' singularities 
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This paper is committed to calculations near a type of future singularity driven by phantom 
energy. At the singularities considered, the scale factor remains finite but its derivative diverges. 
The general behavior of barotropic phantom energy producing this singularity is calculated under the 
assumption that near the singularity such fluid is the dominant contributor. We use the semiclassical 
formula for renormalized stress tensors of conformally invariant fields in conformally flat spacetimes 
and analyze the softening/enhancing of the singularity due to quantum vacuum contributions. This 
dynamical analysis is then compared to results from thermodynamical considerations. In both cases, 
the vacuum states of quantized scalar and spinor fields strengthen the accelerating expansion near 
the singularity whereas the vacuum states of vector fields weaken it. 



It has been determined that our universe is undergoing 
a second era of accelerated expansion [1-5] which can be 
explained by postulating the existence of a fluid whose 
ratio of pressure to energy density must be a negative 
number (see [6, 7] for reviews). This ratio is named the 
equation of state parameter w := pj p and its effective 
value, after including the contributions from all fluids in 
the universe, must be smaller than —1/3 to ensure the 
mentioned accelerated expansion [8]. Fluids with nega- 
tive w are dubbed dark energy in general. The particular 
cases w = — 1 and w < —\ are called cosmological con- 
stant and phantom energy, respectively. 

Phantom energy is disfavored by theoretical reasons: 
it violates all energy conditions [8, 9]. Yet, it is com- 
patible with current observations albeit marginally [10]. 
Because the energy conditions are not satisfied, the sin- 
gularity theorems do not apply and dark energy can lead 
to singularities. Depending on the model, phantom en- 
ergy can produce future singularities that can become 
the end of the universe: the Big Rip singularity [11, 12]. 
This paper deals with softer singularities: the so-called 
sudden singularities [13]. This singularities are softer in 
the sense that divergence occurs in the second derivative 
of the scale factor (see Table I). Moreover, finite objects 
passing through it are not crushed by tidal forces [14]. 

Some effort has gone into investigating whether the 
contributions from the vacuum states of quantum fields 
would enhance or weaken such phantom-energy singular- 
ities [15-20]. A word of caution is needed here. While it 
is possible to exploit the trace anomaly to compute the 
vacuum contributions (at least for conformally invariant 
fields in conformally flat spacetimes) [21], it is not suit- 
able to talk about avoiding the singularity altogether. 
This would imply feeding these contributions, which are 
of first order in ft, back to Einstein's equations. This 
cannot be done consistently in the models we are study- 
ing because we know only the classical behavior of dark 
energy, i.e. up to zero-th order in ft. 

The quantum fields considered in this paper are con- 
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formally invariant and their vacuum state is obtained by 
conformally transforming the Minkowski vacuum. For 
brevity, we will refer to them simply as "fields" but the 
mentioned adjectives are implied. On the other hand, we 
will use "fluids" to refer to classical components. 

The flrst section of this paper is dedicated to flat 
FLRW with non-interacting barotropic fluids. Emphasis 
is given to some monotonicity results applicable if there 
is one dominant fluid. The second section links the be- 
havior of the scale factor to the behavior of the equation 
of state parameter of a dominant barotropic fluid near 
a sudden singularity. It is proven that such link is nec- 
essary and sufficient. In the third section, we compute 
the effects of semiclassical vacuum states of conformally 
invariant fields. There, we show that the quantum en- 
hancing/softening depends on the spin of the field. In 
the last section, general thermodynamic arguments show 
that the heat exchanged by these fields cannot be ne- 
glected. Thus, the models studied in this paper provide 
an example where the thermodynamics of gravitationally 
interacting systems cannot be ignored. We find that the 
dynamical and thermodynamical analysis coincide. Fi- 
nally in the appendix, while sketching a calculation from 
section II, we visit some calculations of a model proposed 
in the literature [16]. 



I. FLAT FLRW WITH BAROTROPIC FLUIDS 

The spatially flat FLRW metric (with scale factor a) 
ds2 = - dt^ + a^(t) ( Ax^ + d?;2 + dz^) 
yields the following Einstein equations: 

2ij--«;2^(p,+K) , (lb) 

where H :— - is the Hubble rate, pi and pi are the en- 
ergy density and the pressure of the z-th fluid, and the 
gravitational constant is set to 1 (k := \/8n). If the 
fluids can be modeled as barotropic perfect fluids, then. 
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by definition, their pressures are functions of their corre- 
sponding energy densities. Let us denote the negative of 
the enthalpy by /, so that 



fi{pi) ■■= -{Pi+Pi) 



(2) 



are functions of p, alone. In terms of the equation of 
state parameters can be written as 



Note that the i-th fluid corresponds to phantom energy 
as long as fi is positive definite. Using Eq. (2), we can 
remove the pressures from Einstein's equations and from 
the conservation of energy equations of each fiuid, which 
now read 



3i7 



Pi + Qi 



fi 



(3a) 
(3b) 



The Qi account for the transfer of energy between fluids. 
As usual, conservation of energy and Einstein's equations 
are not independent. 

Now, wc will make the assumption that there is no 
interaction amongst the fluids, i.e. Qi = 0. This step 
will be justified by its compatibility with the results: If 
at the end of the calculations there is only one dominant 
fiuid, phantom energy near the singularity, whose energy 
density is much larger than the energy densities of all 
other fiuids combined, then we can neglect the interaction 
between fluids. In this case, we can use Eq. (3a) to find 
Pi as functions of a, and Eq. (la) to find the evolution of 
a in time: 



31og(^ 



^- dpi 
da 



(4) 

(5) 



The subindex zero indicates the value of the quantity at 
a given time, say today. 

From Eq. (4), we can reason that pi are monotonic 
fimctions of a as long as fi don't cross zero, cross infinity 
or jump branches {fi could be multivalued). Analyzing 
Eq. (5), we can deduct that a is a monotonic function 
of t provided that the total energy density ^ pi docs not 
vanish or diverge. We can conclude that pi are monotonic 
functions of time until /, or ^ pi vanish or diverge or fi 
change branches. 

Finite-time future singularities will appear if for some 
reason the integral in Eq. (5) cannot keep increasing; so 
that time, in the left hand side of equation of Eq. (5), 
wouldn't keep increasing either. A classification of finite- 
time future singularities was introduced by Nojiri et al. 
[16] and completed by Copeland et al. [7]. It is summa- 
rized in Table I. 

If the integrand in Eq. (5) cannot be evaluated, say 
because the integrand fails to be a real number for a = Ug, 



TABLE I: Classification of finite-time future singularities 
[7, 16]. A dash indicates non-specified behavior. The equiva- 
lences of behavior, p ~ |d| and \p\ ~ |a|, follow from Eqs. (1). 
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then we might have a sudden or type-Ill singularity at 

a,,. On the other hand, if a, the upper limit of Eq. (5), 
can go to infinity, then we have a Big Rip singularity if 
the integral converges; and there is no singularity if the 
integral diverges. 

One corollary out of these statements is that a singu- 
larity type IV cannot be produced by a single barotropic 
non-interacting fluid in a realistic model. This case, al- 
though mathematically possible, is physically impossible. 
Baryonic matter will always be present and its contribu- 
tion to the right hand side of Eq. (la), which evolves as 
o~^, does not vanish because Ug remains finite per defini- 
tion of type- IV singularity. If the energy density of phan- 
tom fluid vanishes then the energy density of baryonic 
matter takes over and drives the evolution of a (without 
any singularity whatsoever). 



II. SUDDEN SINGULARITIES 

A sudden singularity occurs when the pressure diverges 
but the energy density and the scale factor remain flnite. 
By Eq. (1), the first derivative of the scale factor also 

remain finite but the second derivative diverges. Hence, 
a can be approximated, near the singularity, by 



ait) 



,-A{ts-t)+B{tg-ty+^+0{{tg-t) 



(6) 

with positive definite both 5 and e. Possible higher order 
terms in ts — t have been omitted (e > 0) because we 
only need to show the divergence of the second derivative 
while keeping the first derivative finite ((5 > 0). It is easy 
to prove that the behaviors of the total energy density 
and total pressure near the singularity are given by 



12 Pi 



3^2 6AB{6 + 2){tg-t)^ 



E 



Pi ~ - 



2B{d 



a/ (6 + 1) k2 



and (7a) 

(7b) 



Note that the signs in front of A and B in Eq. (6) have 
been chosen in such a way that if A and B are both pos- 
itive then a approaches from below, p also approaches 
Ps from below, and p diverges to — oo. If A was negative, 
then a'{ts) would be negative. But this cannot happen 



3 



because of the monotonicity of a -Einstein equations for 
flat FLRW don't have curvature terms to change the sign 
of a'{t) and we know that a'{t) is positive today. 

Now, let us analyze the conditions under which the 
contributions from fluids like dark matter or electromag- 
netic radiation would not be significant near the singular- 
ity. We need the total energy density to be much larger 
than the dark-matter energy density. This is, the con- 
stant term in the behavior of ^ , cf. Eq. (7a) , must be 
much larger than p^^a^^ / a^^ , where is the energy 
density of dark matter measured today. Therefore, 
and A must be big enough to satisfy 



a2 ^ 3 

A ds ^ (Iq pmo ■ 



If A vanished, then this inequality could not be satisfied. 
This is, only positive definite A is compatible with the 
assumption of non-interacting fiuids. 

If B was negative, the fliiid causing the singularity 
would have positive both energy density and pressure. 
Thus, such fluid would, near the singularity, simply not 



be phantom energy. 

Assuming that only one fluid, dark energy, contributes 
significantly near the singularity, then Eq. (2) for such 
fluid has the form 



/ « sign (B) 



(3 \A\y{2 \B\ {6 + 2))i+* 
K^s + 1)2+5 a 



1+25 



3^2 



(8) 

Although it has been argued that only positive A and B 

are physically relevant for phantom-energy driven future 
singularities, the above formula shows the correct sign 
dependence should A or B he negative. Keeping the 
sign dependence allows for easier comparison with other 
calculations presented in the literature. 

Conversely, if phantom energy is modeled by 



/ = 



C 



+ 0{{p,-py-') 



[Ps-Pf 



(9) 



with positive C, ps, and 5, then the evolution is such that 
the scale factor near the singularity is given by 



a{t) 



1-T + 



3C(l+5) 

— — 



+ 



2{2 + 5)ps 



(10) 



where 



and 




= min 1 2, 



(5-1-3 
J+l 



- 1 - 



1 + ^ 



(11) 



(12) 



Note that Eq. (10) is of the form of Eq. (6) and that 
€ is positive for 5 > 0. While Copeland et al. showed 
that the first term of Eq. (9) yields a sudden singular- 
ity (see Eq. (461) in [7]), the calculation shown here is 
more general in that it only analyzes the behavior near 
the singularity (hence the operator O and the need to 
keep track of e). Thus, it encompasses other models, e.g. 
model (32) in [16], that might behave differently far from 
the singularity. As implied in remarks by Cattoen and 
Visser [22], proving that Eq. (8) can be obtained from 
Eq. (6) is relatively trivial. However, the proof of the 
converse is rather laborious because one must prove that 
the dismissed terms in Eq. (9) can also be dismissed in 
Eq. (10). Such proof follows the lines of the calculation 
in the appendix. 

We reach then the following conclusion: a phantom- 
energy model where barotropic dark energy is the only 



TABLE II: Spin-dependent coefficients in Eq. (14). 
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significant fluid near the singularity will produce a sud- 
den singularity, Eq. (6), if and only if its behavior near 
the singularity has the form of Eq. (9). The relationship 
between {A,B) and (C,ps) can be read off from equa- 
tions (6) and (10). One implication of this theorem is 
that sudden singularities cannot be achieved with a static 
equation of state parameter, it must evolve according to 



C 



w 



P{Ps 

near the singularity. 



pY 



= 0(r-'+T+7) 



(13) 



III. SEMICLASSICAL FIELDS 

We will use the semiclassical expression for the vacuum 
stress-energy of conformally invariant quantized fields in 
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a vacuum state conformally obtained from Minkowski 
spacetime [21]: 

{Tiiv) = (^iiv R'^-a ~ R,fif + RR^v — '^9lJ.v R^ + 
/3 i^RRnv — RJi Rva + 2^nv Rar R'^^ — -^9nv R^ , 

(14) 



where R is the Ricci scalar, i?^^ is the Ricci tensor, and 
a and (3 are spin-dependent coefficients given in Table II. 
Because the derivatives of R diverge faster than R or g, 
the first two terms of the coefficient of a are the ones 
that contribute the most: 



4 

K p. 



Pa ■■= (roo)|spi„=a = a^{3C)Tt^6{{l + 6)T)-^+^ , 
Pa ■■= (Tn)|spin=a = ""^(^ C)rh S {1 + 2 S) a,^ ((1 + S) t)-'+^ . 



(15a) 
(15b) 



r 



Both expectation values diverge as the singularity ap- 
proaches, the pressure faster than the energy density. 
The equation of state parameter of the quantum con- 
tributions is then given by 



w ■ 



, 1 + 2S 

3(1 + 6)' 



(16) 



It is negative, it doesn't depend on the spin of the field, 
and it diverges at the singularity faster than Eq. (13). 
These divergences mean that the approximation breaks 
down before the singularity occurs. Nevertheless, a qual- 
itative analysis is in order. The behavior of the pertur- 
bation is determined by the sign of a, which is the only 
quantity not defined positive in Eq. (15b). Because a is 
positive for both scalar and spinor fields, then pa > 0, 
Pa < and therefore these vacuum states enhance the 
singularity. This happens because adding these vacuum 
perturbations is equivalent to adding more phantom en- 
ergy to the right hand side of Eq. (1). Vector fields, with 
negative a, counter the contributions from dark energy 
and therefore soften the singularity. 



IV. THERMODYNAMICAL CONSIDERATIONS 

While the vacuum state of a system described by an 
exact Lagrangian has no thermal properties, it is not 
impossible for a vacuum state to become thermal (e.g. 
[23, 24]). In general (see [25] for a canonical exposition), 
small unknown terms in the Hamiltonian cause the ther- 
malization of the system (equation (12.2) leads to equa- 
tion (13.13) in the reference). 

In this section, the vacua are considered as individ- 
ual subsystems. We admit that their evolution is known 
only up to some terms of first order in h and that their 
complete Lagrangian would contain higher order terms 
in h, self- interactions, and interactions with the phan- 
tom fluid. Thus, these vacua qualify for thermal calcu- 



lations. The previous section provide the energy density 
and pressure of the fields. The pressure is interpreted as 
the partial pressure of the vacuum state in the mixture 
of cosmological fluids. We don't need to assume that the 
process is quasi-static because we are not concerned with 
the temperature or the entropy of the subsystems; we 
shall be satisfied with computing the sign of the change 
of the enthalpy of formation. 

In the first law of thermodynamics 6U = SQ — 5W, we 
can replace 



U (X paa^ ; 

W (XPatt^ ■ 



(17a) 
(17b) 



The approximations for Eq. (14) are valid only for free 
fields (see [26] page 5). Thus, the lack of interaction 
terms in Eq. (17b) is justified. We are interested in the 
sign oi SQ because it determines whether the expansion 
is exothermic or cndothermic. As shown in Eqs. (15), 
the pressure diverges faster than the energy density and 
therefore 6Q « 6W. Hence, 



sign {6Q) = sign (3pa Sa + Spa) 



(18) 



We can now dismiss the first term in the right hand side 
because it behaves as r i+'S which is slower than the 
T~'*^T+^ divergence of Spa- Therefore, the sign of the 
heat flowing into the system is the same as the sign of 
pressure change. The latter is determined by the neg- 
ative of the sign of a because all of the other quanti- 
ties in Eq. (15b) are positive definite. The expansion 
is then exothermic for both scalar and spinor fields and 
endothermic for vector fields. Exothermic reactions are 
spontaneous and thus they enhance the singularity. We 
conclude again that scalar and spinor vacua enhance the 
singularity and vector vacua soften it. This is in agree- 
ment with the dynamical results of the previous section. 

While the models considered in this paper furnish an 
example of dynamics and thermodynamics coinciding in 
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their predictions, this might not necessarily be true in 
general. A counterexample would be extremely valuable 
because it will point to a deficiency of Quantum Field 
Theory in Curved Spacetimes (QFTCS). This theory is 
considered complete up to first order in H. But if the 
dynamics and thermodynamics don't coincide, then the 
predictive power of QFTCS would be limited because it 
wouldn't be possible to determine the evolution of an 
entropy-driven system without a Statistical Mechanics 
theory of gravity (let alone Quantum Gravity). 



APPENDIX: ANALYSIS OF MODEL (32) OF [16] 

In this appendix we will follow the steps necessary to 
obtain Eq. (6) from Eq. (9). The general case is rather 
complicated because of the need to keep correct track of 
the orders of magnitude, and at some point it involves ex- 
panding a hypergeometric function composed with a log- 
arithm evaluated at the singularity of the logarithm. So, 
instead of cluttering with long mathematical expressions, 
we will perform the basic steps while reviewing calcula- 
tions in a model proposed by Nojiri et al. (model (32) of 
[16]). Also, this derivation will expose some problematic 
points in case the reader is interested in reproducing the 
full calculations. 

Consider dark energy modeled as a fluid with 



fiP) 



(A.1) 



where Ps > is the dark energy density at the sin- 
gularity [27], and 7^0. If ps was negative, then 
/ might not be defined over the real numbers. Using 
the monotonicity results from section I, it is trivial to 
show that Eq. (A.l) corresponds to phantom energy as 
long as and whenever po < ps- On the other hand, if 
po > Ps, the fluid described by Eq. (A.l) might begin its 
evolution as non-phantom dark energy (for 7 > and 

Po > + ^1 + Ah/^p^'^ ) /2) or a normal fluid (i.e. 
positive both pressure and energy density) but it will be- 
come a normal fluid near the singularity. 
With this model, Eq. (4) yields 



Po') 



(w) 



(A.2) 



Even though the above equation might have two signs at 
front of the radical sign, one of them is spurious. The 
sign shown yields p ^ po when a ^ uq- 

As argued previously, the singularity occurs when the 
quantity inside the square root vanishes. This defines the 
value of a at the singularity: 



uq exp 




(A.3) 



If b is negative, then the radius at the singularity will be 
smaller than the radius today. In that case, the contri- 
bution from matter and dark matter is not diluted and a 
model without such contribution is unphysical. We will 
continue our analysis assuming 6 > 0. 

Since we are interested in the behavior near the singu- 
larity, instead of Eq. (5), Eq. (la) can be approximately 
written as 



%/3 



{ts - t) 



da 



iVpia) 



(A.4) 



The contributions from other fluids but dark energy have 
been discarded on the basis that as is much larger than 
oq. The integrand can be expanded around Us, the in- 
tegral can be evaluated, and the resulting series can be 
inverted. The first three terms of the outcome are 



a{t) = Us [ 1 — T 



Ti]+0 (r^) (A.5) 



where A :— sign (log(po/Pis)) and r as in Eq. (11). 

Note that the second derivative of a is the first to di- 
verge. Thus, the singularity is only a sudden singularity, 
no matter the value of 7. Contrast this to the rich 
7-dependent structure reported elsewhere [16]. 

The independence from 7 could also be deduced from 
the behavior of Eq. (A.l) near the singularity: 



bp 



2(1-7) 



l^\Ps- P\ 



+ 0{{ps-pf) (A.6) 



which shows that ps is a single pole and therefore falls 
within the theorem of section II. 
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